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Abstract. Tropical R is the birational map that intertwines products of geo- 
metric crystals and satisfies the Yang-Baxter equation. We show that the D„ 
tropical R introduced by the authors and its reduction to and are 

equivalent to a system of bilinear difference equations of Hirota type. Asso- 
ciated tropical vertex models admit solutions in terms of tau functions of the 
BKP and DKP hierarchies. 



1. Introduction 

Let B — {x — (x\, . . . ,x n )} be the set of variables. The tropical R for A 



(i) 



[HHIKTT, Ki, Y] is the birational map R : B x B 



(l.l) 



R(x,y) = (x ,y), 

n n 

P l (x,y) = Y / l[ x 
fc=l j—k 



B x B specified by 

P l - 1 (x,y) 



y, 



Pi(x,y) 



' -J 1 I Vi+j: 
3 = 1 



where all the indices are considered to be in Z/nZ. It satisfies the inversion relation 
R 2 — id on B x B and the Yang-Baxter equation R1R2R1 = R2R1R2 on B x B x B, 
where R\(x,y,z) — (R(x,y),z) and R2(x,y,z) = {x,R(y,z)). Given (x,y), it is 
characterized as the unique solution to the equations on (x' , y'): 



(1.2) 



i iji 



1 



1 



1 
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with an extra constraint T[™ =1 {xi / yl) = YYi^iiVi/Xi) = 1- 

A representation theoretical background for the tropical R is provided by the 
geometric crystals [BK] and their natural extrapolation into the affine setting as 
demonstrated in section 1 in [KOTY]. To recall it consider the matrix realization, 
where each element x € B is associated with the matrix (see also [KNY]) 



M(x,() = 



-1 



V 



-1 x- 1 / 



involving the spectral parameter ^. Then the structure of geometric crys- 

tal on B is realized as simple matrix operations. For instance the geometric 
Kashiwara operator el : {x\, . . . , x n ) 1— »•(... , Xi-i, cXi, c~ 1 Xi+i, Xi+2, ■ ■ ■) is induced 
by a multiplication of M with certain unipotent matrices. The product of ma- 
trices M(x^\ C) ' ' ' M(x^\ C) corresponds to the product of geometric crystals 



(xW, . . . ,xV)) e B x • • • x B in the sense of [BK]. Eq.(1.2) is equivalent to the 
matrix equation 

(1.3) M(x, ()M(y, C) = M(x', QM(y', (). 

Due to the presence of the spectral parameter £, its non-trivial solution is unique, 
which characterizes the tropical R as the intertwiner of the products of geometric 
crystals. 

Geometric crystals are so designed that the relevant rational functions become 
totally positive [BK, BFZ, L]. There is no minus sign in Pi(x,y) in (1.1) indeed. 
Such functions can consistently be transformed into piecewise linear ones by replac- 
ing +, x and / by max(min), + and — , respectively. Consequently, the structure of 
geometric crystals reduces to the one for crystal bases [K]. In the present case, the 
expression (1.1) leads to the piecewise linear formula [HHIKTT] for the combinato- 
rial R [NY]. An analogous result is available in [KOTY] for the combinatorial 
R [HKOT]. The above transformation can be realized as a certain limiting pro- 
cedure [TTMS] and is often called the ultradiscretization (UD). As related topics, 
we remark the tropical combinatorics [Ki, NoY] and the soliton cellular automata 
associated with crystal bases [TS, TNS, HKT1, FOY, HHIKTT, HKOTY, HKT2] . 

Now let us turn to the aspect of the tropical R as a classical integrable system, 
which is the main subject of this paper. Recall the discrete time Toda equation 
[HT, HTI] on the electric "current" I\ and "voltage" V*: 

with periodicity If — J* +n , V* — V* +n . This difference equation is known to be 
integrable. Eq.(1.2) is identified with it via x~ x — J| +1 , y^ 1 = V*, x'J 1 = V t t+1 and 
y'7 t = [Y]. The matrix equation (1.3) is a Lax representation in this context. 
In fact the tropical R is equivalent to a system of bilinear difference equations of 
Hirota type, which was effectively the base of the analyses in [HHIKTT]. To see 
this, introduce the functions rf (1 < J < 4,i e Z/nZ) and the parameters \,Ki, 
and make the change of variables 

x^ 1 = XiSrf/Srf, y^ 1 = k^t^/StI, 

with St/ — t/ /t/_ 1 . Then the former relation in (1.2) is automatically satisfied 
and the latter is translated into 

(1.4) \rtx< - Kirfrl, = arlrl, 

for any nonzero parameter a independent of i. The birational map R : (x, y) i— > 
(x',y') is induced by an automorphism rf <-» rf, A, <-> Ki, a — > —a of (1.4). Eq. 
(1.4) is a version of so-called Hirota-Miwa equation [H, M], a prototype discrete 
soliton equation to which the well-developed machinery of free fcrmions and infinite 
dimensional Lie algebras [DJM, JM] can be applied. The solutions are provided 
by tau functions of the KP hierarchy with a certain reduction and discretization of 
time variables. 

In this paper we elucidate the classical integrability of such sort for the DIP 
tropical R introduced in [KOTY] together with its reductions to ^n-i an d Cn^ ■ 
We show that they are equivalent to a quartet of A type bilinear equations like 
(1.4), and construct solutions in terms of tau functions in soliton theory [JM]. The 
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cases Dn and A^-i ar e related to reductions of the DKP hierarchy associated to 
the algebras and D 1 ^ . A key is the relation between two kinds of tau functions 
that originate in the isomorphism ~ D'^. See Lemmas B.l and B.2. The 
case corresponds to a further reduction to the BKP hierarchy associated with the 
algebras ~ B'^ . 

We shall formulate tropical vertex models, where the tropical R plays a role 
of local time evolution. Namely it is the two dimensional system on a square 
lattice where each edge is assigned an element of a geometric crystal in such a way 
that those four surrounding a vertex are related by the tropical R. Our solutions 
to the bilinear equations are naturally extended to the tropical vertex models by 
duplicating the original lattice and attaching a Clifford group element to each face 
of it under a certain rule. The Yang-Baxter equation for the tropical R is naturally 
understood from such a point of view. See Remark 5.6. 

We expect that the results in this paper are glimpses of deeper relations yet to be 
explored between geometric crystals, crystal bases, solvable lattice models, discrete 
and ultradiscrete solitons and so forth. A rough picture at present looks as follows. 



[/^-modules 
Quantum R 
Vertex models 



Crystals 
Combinatorial R 
Soliton cellular automata 



UD 



Geometric crystals 
Tropical R 
Soliton equations 



In [D], Drinfcld proposed set-theoretical solutions of the Yang-Baxter equation 
as one of unsolved problems in quantum group theory. Since then a number of 
approaches have been made. See for example [ABS, ESS, JMY, MV, O, VI, V2, 
WX] and references therein. Our tropical R is a distinguishable example which 
is placed in the unique spot in the above picture and enjoys the total positivity. 
It should be noted that the combinatorial R is also a remarkable example of the 
set-theoretical solution which intertwines products of finite sets. 

The paper is organized as follows. In Section 2 we recall the tropical R for 
following [KOTY]. In Section 3 the bilinear equations are introduced, which are 
divided into four groups corresponding to each vertex in Figure 1. Uniqueness and 
existence of the solution are shown in Proposition 3.5. In Section 4 we establish 
the bilincarization of the tropical R in Theorem 4.4. In Section 5 we construct 
solutions of the bilinear equations in terms of tau functions of the DKP hierarchy 
formulated with two component fcrmions. The final result is given in Theorem 
5.5. In Section 6 we introduce the tropical R for -A^n-i an d Cn^ as reductions of 
the D ( n ] case. Parallel results on the bilincarization and solutions arc obtained. 
Appendix A provides a proof of Lemma 4.3. Appendix B contains elements of the 
free fcrmion approach [JM] as well as the key lemmas B.l and B.2. We leave the 
calculation of the ultradiscrete limit of tau functions as in [HHIKTT] for A^li as 
a future problem. 

2. Tropical R for £>i 1} 

Let us recall the tropical R for from [KOTY] . Let x = (xi, . . . , x n ,x n -i, . . . , a?i) 
be a set of variables. The geometric £>„^-crystal is the set B = {x} equipped with 
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additional structures such as tpi, Si and the geometric Kashiwara operators e\. The 
product B x • • • x B again becomes a geometric crystal, which is an analogue of 
the tensor product of crystals. The tropical R is a birational map B x B — > B x B 
commuting with the geometric Kashiwara operators. Leaving the precise definition 
to [KOTY], we concentrate here on its explicit form. 
Set 

(2.1) £(x) = X\X 2 ■ --XnTn-i ■ ■ -X 2 Xi, 

and call it the level of x. On B we introduce the involutions o\ and er„ which preserve 
the level. Explicitly for x = (x\, . . . , x n ,x n -\, . . . , xi) <G B, they are defined by 

(2.2) (7 1 : xi < — > x\, 

(2.3) (7 n . X n —\ ► X n ^\X ni X n —l > X n —\X ni x n > \jx n ^ 

where the variables not appearing in the above are unchanged. On B x B we 
introduce o\ , a n analogously as 

(2.4) cr a (x,y) = (<7 a (x),cr a (y)) a=l,n. 
Furthermore we introduce the involution by 

(2.5) a*:xi< — >y i7 Xi < — > yi (l<i<n-l), x n < — > y n 

for the elements x = (x\, . . . ,x\) and y = (y\,...,y 1 ) of B. Note that ai,a n 
and o * are all commutative, among which c* is the only one that mixes x and y 
and interchanges the levels. Given a function F = F(x,y), we will write F° a = 
F(a a (x,y)) for a=l,n, *. 
We set 



n — 1 /m— 1 \ / \ n—1 

(2.6) Vo =i(x)^+i(x)j2 n - i+^)+^n 

n—1 /m— 1 \ / \ n—1 

+<+w£ (n|)( 1+ £) + n^, 

and define Vi (1 < i < n — 1) by the recursion relations: 



(2.7) 



Vi = y> (=Vi-x + (£(x) - i{y)) (- + ±S) , (1 < i < n - 2) 
K-r = ^ (^-K- 2 + (/(*) - + ^-)) • 



In terms of Vi we also define Ui (1 < i < n — 1) by 

















2/i 



The tropical R is a birational map R(x, y) = (x 1 , y') on B x B given by 



Vq 1 _, _ V 



(2.9) 



, Vi-iUi _, _ Vi_i 
^ = ^7777 — » ^ = 2/;——, (2<*<n-l) 

I _ Vn-1 
x n — Vn -, /CT „ i 



CTl 



yd, y<J, JJ_ 

y'i = XiT0, y\ = xi 1 , (2<i<n-l) 



yvn 
v n—1 

Vn-1 



In [KOTY], C/i here was denoted by Wi in Definition 4.9, and their transformation 
property under a a was summarized in Table 1. As shown therein, the tropical R 
is subtraction-free (totally positive), interchanges the levels, commutes with the 
geometric Kashiwara operators, and satisfies the inversion relation R(R(x,y)) = 
(x, y) and the Yang-Baxter relation. 

Remark 2.1. Let (x',y') = R(x,y) under the DrP tropical R. It is easily seen 
that x n = y n — 1 is equivalent to x' n — y' n = 1. Similarly, (x\,yi) — (xi,yi) is 
equivalent to (x' 1 ,y' 1 ) = (x^yi). 



3. Bilinear equations 

Here we introduce a system of bilinear equations and study its properties such 
as existence and uniqueness of the solution. It is a preparation for Section 4, where 
the equations are related to our tropical R. 

In the rest of the paper we fix the elements A = (Ai, . . . , Ai) and k = (k\, . . . ,K\) 
of B, and set 

(3.1) l = t(X), k = £{n). 

We assume lk(l — k) ^ throughout. We introduce the variables which we call the 
tau functions: 

Si,Wi,Ni,Ei (l<i<n-2), t?,t/ ,rf ,rf ,rf (0 < * < n). 

They will be called generic if they are nonzero. Let a and be arbitrary nonzero 
parameters. Our bilinear equations read as follows: 

MiViT 4 - Al^iTo 2 = ot?tI, 
KiE^Ni - \ l N t _ 1 E l = ar^rl (2 < i < n - 2), 

K n -iK n E n -2Tn - A„_iA„Y„_ 2 T I t = aT^T*, 
K n -lE n - 2 Tn-l - K-l N n - 2 T^_ 1 = (XT^T^^. 
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The i-th equation here will be referred as (0 < i < n). 

KxN^f + \ 1 W 1 tI =/3r °r 2 , 
^i-ZViTq + AiM^iTq = /3t^t 2 , 
KiWi-iNi + XiNi-iWi = 0t?t? (2 < i < n - 2), 

KnKn-lWn^Tn-x + X n -1 A r „-2T^_ 1 = /?T^_ 1 T^_j, 
Kn-lW n -2T^ + A„_iA„iV„_2T^ = /3r°r 2 . 
The i-th equation here will be referred as (2, i) (0 < i < n). 

KlSn$-\lWlT$=OT?Tf, 

KiW^Si - XiSi-tWi = onfif (2 < % < n - 2), 

The i-th equation here will be referred as (3, i) (0 < i < n). 

KiSiTi 1 + Ai^ri 3 - /3r °r 4 , 

Ki5ir 1 + Ai£;i7i? = /9TM. 
KiEt^Si + XiSi-xEi = fa^rf (2 < i < n - 2), 

K n -lE n -2Tn-l + A n A„_iS'„_2T^_ 1 = /3t°_ 1 T*_ 1 , 

k„-ik„-E„_2T^ + A„_iS„_2T* = /3r°r^. 

The i-th equation here will be referred as (4, i) (0 < i < n). 

We call {(J,i) | 1 < J < 4, < i < n} the Z)^ bilinear equations. 

Lemma 3.1. For generic input Ni, . . . , N n -2, W\, . . . , W n -2 and Tq , . . . ,t£ with 
J = 1, 2, 3, there exists a unique solution to the equations ( J, < J < 4, 1 < i < 
n — 1) on f/ie variables Si, ... , S n -2, Ei, ... , E n _ 2 , t®, . . . , t°_ x and Tq , . . . , t*. 

Proof. The equations (2, 1) — (2,n— 1) fix r^, ... , t^_ 1 . Regard and (3, i) with 
l<i<n— las linear equations on the unknowns Si, ... , S„_2, E\, . . . , E n - 2 , Tq 
and t*. Arrange them in a matrix form Ax = ah, where the column vectors x and 
b are specified by 

x = *( T > Ei,.. ., E n - 2 ,T*, S n -2, ■ ■ ■ , Si), 

h _t ( 1 1 1 3 3 3n 

Then one finds that the 2n — 2 by 2n — 2 coefficient matrix A is an upper triangular 
Jacobi matrix except the bottom left clement — AiWi. Hence it is easy to see 
dct A = (k — 1)tqT% nr=i which is nonzero for generic input. Therefore the 

unknowns in x are uniquely determined. With rf , . . . , r°_ 1 and x at hand, one can 
determine the remaining ones rf , . . . , t^_ x from (4, 1) — (4, n — 1). □ 

According to Lemma 3.1, the bilinear equations ( J, 0) and (J, n) are not needed 
to determine the tau functions that will appear in the parameterization (4.2)-(4.3). 
However, they are essential in order to incorporate the involutions a a in the previous 
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section into tau functions. The rest of this section concerns this point, leading to 
Propositions 3.5, 3.6, and ultimately to Proposition 4.2. 

Lemma 3.2. Suppose that St, W t , N h E t (1 < i < n - 2) and t? (0 < J < 4, < 
i < n) satisfy ( J, < J < 4, 1 < i < n — 1). Then we have 



(3.2) 

(3.3) oltI_ x tI + pXnT^T^ - a\ n K n rl_ x rl 

Proof. Write (1,1) — (4, 1) in the matrix form: 







V^ITq 3 



Ait 3 









KlT 2 






_0 

-Air 4 



\ 






/orA 




Ei 




/3r 4 




Si 




3 


/ 


\ W1 J 







The matrix on the left hand side is annihilated by multiplying the row vector 
(AiKiTq , AiKiTq , — AiKiT 1 , — AiKiTq) from the left. So the same should happen on 
the right hand side, proving (3.2). The relation (3.3) is verified similarly by using 
(l,n- 1) - (4,n- 1). □ 

Lemma 3.3. Suppose that Si, Wi,N h E, t (1 < i < n-2), rf (1 < J < 4, < i < n) 
and rf (1 < i < n — 1) satisfy the equations ( J, < J < 4, 1 < i < n — 1). Then 
there exists a unique solution to the equations (J, 0) (1 < J < 4) on Tq . 

Proof. It suffices to check compatibility of the four equations ( J, 0) (1 < J < 4), 
that is, they all lead to the same Tg. Multiply (1,0) by Tg, and (3,0) by r 1 . The 
difference of the resulting left hand sides can be grouped as 

Due to (2, 1) and (4, 1) this vanishes, proving the compatibility of (1,0) and (3,0). 
The compatibility of (2,0) and (4,0) is similarly shown by using (1,1) and (3,1). 
Let us check the compatibility of (2,0) and (3,0). We are to show 

ar^^N^l + \ 1 W 1 t1) - fSr^mS^ - A^r 4 ) = 

Upon multiplying Ki, the first term contains the factor KiN^ and the third does 
kiSiTq . Eliminate them by using (2, 1) and (3, 1), respectively. After a cancellation 
the resulting expression becomes W\ times the left hand side of (3.2). □ 

Lemma 3.4. Suppose that Si, W u N u E t (1 < i < n-2), rf (1 < J < 4, < i < n) 
and Tj (1 < i < n — 1) satisfy the equations (J, < J < 4, 1 < i < n — 1). Then 
there exists a unique solution to the equations (J, n) (1 < J < 4) on t°. 

Proof. This is verified similarly to Lemma 3.3 by using (3.3). □ 

Since the compatibility of the two systems (1,0) — (4,0) and (l,n) — (4, n) is 
trivial, Lemmas 3.1, 3.3 and 3.4 lead to 

Proposition 3.5 (Unique existence). For generic input N\,..., iV„_ 2 , W\, . . . , W n -i 
and Tq , . . . , t£ with J = 1,2, 3, there exists a unique solution to the bilinear equa- 
tions (J,i)(l < J < 4, < i < n) on the remaining variables S\,...,S n -2, 
E u ... , E n _ 2 and t j , . . . , r r { (J = 0, 4). 
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cri : 


r '«— 


ri 7 (0<J<4), 


0~n : 


r n-l ^ 


-r„ J (0<J<4), 


(7* 


■ r n-l ^ 


o 2 : 






r? (0< i <n-2), 



Let us proceed to the automorphism of the bilinear equations. For brevity the 
array T = (A, k; Si, W%, Ni,Ei,rf) consisting of all the tau functions together with 
(A, k) G B x B will be called data. We extend the involutions a a (a = l,n, *) on 
B x B defined in the previous section to the data as follows: 

(3.4) a a (a = l,n, *) acts on (A, k) e B x B by (2.2)-(2.5) as before, 

(3.5) 

(3.6) 

(3.7) <7* : r°_ 1 < ► t°, t^_ x < ► r^, r^_j < ► 

3 l 3 

VF.^iV,, Si^Ei (l<*<n-2). 

In each involution, the tau functions not appearing in the above are unchanged. 
See also Proposition 4.2. It is immediate to check 

Proposition 3.6 (Invariance). The involutions a a (a = l,n, *) defined in (3.4)- 
(3.7) on the data T are automorphism of the bilinear equations (J,i)(l < J < 
4,0 < i < n). 

In particular, the composition C7„<7* interchanges and (3,i) leaving (2,i) 

and (4, i) invariant. 

4. BlLINEARIZATION OF TROPICAL R 

Given <t = (<*i, . . . ,7ii) £ B and the arrays 

C = (Ci, . . . , C„_ 2 ), t = (t , . . . , r„), t' = (tq, . . . ,T' n ), 

we define 

\fJL\T, C,t'] = (zi,...,Z n ,Z n -l,...,Zl) e B, 
-l^OTl _ iTITq 

-1 ClT 2 _ 1 C\T 2 



C2Tqt{ ' C 2 T ri' 
^ = ^9^. ^=W 1 ^ (3<^<n-2), 

_1 C n -2Tn-l _ 1 C n -2T n 

Z n -1 — / i z n-l — Mn— 1 / J 

T„-lT n _ 2 T„_2T„ 

z n — Mn 7 

T « T n-l 

in case n > 4. For n = 3 we slightly modify it into 

-1 T 0T! _ xTiTq 

z i = Mi "Tt 7 z i = Mi ' 

-l C\t' 2 _ 1 C\t 3 

z 2 - A*2 - — TT> Z2 = ^2 7i 



= fJ- 3 



T2T Q T 1 T TiT 3 
T3T2 



Note that £([fi; r, C, r']) = ^(m) -1 - From the data T = (\, k; S i ,W l , N u Ei,r- J ), we 
construct x,y,x',y' £ B as 

(4.2) x= [X;t 3 ,W,t 2 ], y = [«; r 2 , JV, r 1 ], 

(4.3) x' = [k;t 3 ,S,t% y' = [A; t 4 ,E, t 1 ]. 
Note from (3.1) that 

(4.4) £(z) *(a/) = £{y) = k~\ 

Given any pairs (x,y) and (A, k) such that £(x) — ^(A) -1 and £(y) = I?(k) _1 , the 
parameterization (4.2) is always possible due to 

Proposition 4.1. For any z,fj,GB such that£(z) = ^(yu) _1 , the equation [fj,; r, C, t'] = 
z with fixed t, (resp. r' ) on the variables (C, r') (resp. (C,t)) admits a one pa- 
rameter family of solutions. 

It is customary to attach a vertex diagram (cross) to the relation like R(x, y) = 
(x',y') for a quantum or combinatorial R, where the four edges correspond to 
x,y,x',y'. Figure 1 is an analogue of such a diagram, which may be of help to 
recognize the pattern (4.2)-(4.3). Elements of B are represented there with double 
lines. Consequently the vertex diagram is separated into several domains where 
various tau functions live. 

Figure 1. Tau functions 



y 





T 2 


N 


T 1 


X 


W 




E 


y' 




T 3 


s 


T 4 



x' 



It is readily seen that the transformation T — ► er a (T) of the data induces the 
change (x,y) — > cr a (x,y) and (x',y') — > a a (x',y'). More generally we have 

Proposition 4.2. Let F = F(x,y) be any function on B x B and 
T = (A, k; Si, Wi, Ni,Ei,rf) be any data obeying all the bilinear equations. Denote 
by F\q- an expression of F resulting from the substitution of (4-2) followed by any 
possible application of the bilinear equations. Then one has F aa \r = o a (F\-r) (a = 
l,n, *), where a a on the left hand side is specified by (2.2)-(2.5), whereas on the right 
hand side by (3.4)-(3.7). The same fact is valid also for (4-3) and F — F(x' ,y'). 



Proof. This is due to Propositions 3.5 and 3.6 
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□ 



Lemma 4.3. Suppose that the data (A, k; Si, Wi, N i7 E i7 rf) satisfies all the bilinear 
equations. Then under the substitution (4-2), the functions Ui,Vi defined in (2.6)- 
(2.8) are expressed as follows: 



2 4 2 4 22244 

v o- V - "~T~3~' v a ~ 1 ~ x i x i x 3 x 3 ' 

T '0 T l T l T T l T T l 

^T-2_4 -,--2 4 2 2 2 4 4 

t 7 T/ "! Inn T/* 7 * I/* 7 ™ . n ~ _ n ~ 1 /Y - n n— 1 n 

Vn-1 - V n -1 - 13 ' "n-1 - "n-1 ~ 1 3 ' ~ _1 T 1 T 3 T 3 ' 

'n'n ' n— 1 n — 1 ' n—1 n ' n—1 n 

where 7 = (fc — /)/?/ (Ika). 

The proof is available in Appendix A. The main result of this section is 

Theorem 4.4 (Bilinearization of tropical i?). Suppose that the data (A, n; Si, Wi, N i7 

Ei,rf ) satisfies all the bilinear equations. Then x,y,x',y' specified in (4-2) and 
(4-3) obey the relation R(x,y) = (x',y'). 

Proof. Substitute (4.2) into the right hand sides of (2.9) and apply Lemma 4.3. 
Then the result agrees with (4.3). □ 

Given (x,y), the tau functions r 1 , t 2 , t 3 , N, W satisfying (4.2) are not unique. 
Theorem 4.4 guarantees that (x',y') in (4.3) is independent of their choice under 
the bilinear equations. From Proposition 4.2 and the remark preceding it, we find 
that R(a a (x,yj) = a a (R(x,y)) in agreement with [KOTY]. We conclude that the 
tropical R acts on the data as the interchange 

A < — ► k, t?^t?, Wi^Si, Ni^Ei 

for all i. Combined with a — > —a, this is another automorphism of our bilinear 
equations that interchanges (2, i) and (4, i) leaving and (3, i) invariant. This 
is a complementary transformation with the o- n o~* mentioned after Proposition 3.6. 



5. Solutions of bilinear equations 

There is a link between our bilinear equations and the soliton theory in terms of 
tau functions and infinite dimensional Lie algebras [JM]. We show in Section 5.1 
that certain vacuum expectation values of the 2 component free fermions provide 
solutions to the bilinear equations. They are parametrized with the elements of 
the algebras D 1 ^ ~ and their reduction. Then in Section 5.2 the solutions for 
the local equations are naturally extended to the tropical vertex model on the two 
dimensional square lattice where the tropical R plays the role of local time evolution 
at each vertex. As for notations and basic facts on the free fermion approach, see 
Appendix B. 

5.1. Tau functions as vacuum expectation values. In this subsection, the 
letters x, y, etc. stand for the array of infinitely many time variables as in Appendix 
B. They should not be confused with elements in B, which have now been effectively 
replaced by the tau functions. Thus we refresh Figure 1 into the following: 
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Figure 2. Increment time arrays 



UK-*) s(K-i) 



eiL- 1 ) 
i{L-i) 



T 2 


N 


T 1 


w 




E 


T 3 


S 


T 4 



Here we have introduced the parameters K,L, which will play a role of level 
in B. Each line is assigned with the array e or e determined from these pa- 
rameters according to (B.2). We assign the time variables x J on the 9 domains 
J = 0, 1, 2, 3, 4, N, S, W, E in Figure 2 so that the following recursion relations are 
satisfied: 

x — x — x x == x x — s(L 
x 2 - x w = x N - x° = x 1 - x E = e{L- 1 ), 
(5J) x N - x 2 - x° - x w = x s - x 3 = iiK- 1 ), 

x 1 - x N = x E - x Q = x 4 - x s = eiK- 1 ). 

These relations determine all the x J, s consistently upon specifying any one of them 
as an initial condition. We fix x 1 to an arbitrary odd element in the sense of 
Appendix B.l. Then from the construction (5.1), all the time variables on the 
corner domains are odd, i.e., 

(5.2) x~J = x J J = 1,2, 3, 4. 

On the other hand the variables in the domains TV, W, S, E satisfy 

x^ = x N + eiK- 1 ) - e{K- x ), 7 s = x s + e{K- x ) - e^ 1 ), 
X™ = x w + e(L- 1 )-£(L- 1 ), 7 E = x E +s(L- l )-s{L- 1 ). 
For i > 1 define further 

x i = x + z i = { x i,l ? x i,2i x i,3' • ■ 
i 

(5.4) Zl =0, ^-^e^ 1 ) (i>2), 

fe=2 

where a 2 ,a3, . . . are nonzero parameters. Note that xf = x J and xf + e(a^ 1 ) = 

x i-l- 

Let g' e e D °° and g e e D °° be the elements that are related as in (B.5). We shall 
use the vacuum expectation values (B.6) exclusively for odd y = (yi,0,j/3,0, . . .) 

li 



with fixed 2/1,2/3, . . .. Thus we simply write them as Fi lt i 2 - t i(x;g) and fi(x;g r ). Our 
solution to the bilinear equations is constructed in two steps. In Step 1, we construct 
solutions in the infinite n limit. In Step 2, we will impose certain reduction condition 
on the elements g, g' to make the equations close for finite n. 

Step 1. We identify the tau functions with the vacuum expectation values as 

(5.5) X 3 = F hl (xf;g) j>l (X = S, W, N, E), 

(5-6) -J=\ f l { fj'\ ^o' 1 ' U = 1.2.3,4), 

[F 1A (xj;g) j>2, 



(5.7) r) 



o_ J(i^ + (-l) J tfb,i;i)(3?;0) 3 =0,1, 



In the bilinear equations (J, i) in Section 3, consider the n infinity limit, where one 
just forgets (J, n — 1) and (J, n) for all J = 1, 2, 3, 4. We call the resulting system 
the D^) bilinear equations. 

Lemma 5.1. 

*i,i(*i)=ToV J =1,2,3,4, 



2fb,i«;fl) 1 _ r i T 2 ±r i T 2 

2iF 0>1;1 (x?;g) }~^i±^ , 

2F , 1 «;< 7 ) \_ 2 3 , 2 3 

2iF , 1;1 «; 9 ) J- TlTl ±T ° T °' 



2^0,1(^1; 3) I _ 3 4 1 3 4 



2iF ,i ; i(^; 9 ) 
2Fo,i(x?;g) 



"l I _ _4 1 , 4 1 

2zF , 1;1 (xf; 5 ) j"^i±W 

Proof. The first equality follows from Lemma B.l and (5.2). The other relations 
are due to Lemma B.2 and (5.3). □ 

Proposition 5.2. The parameterization (5.5)-(5.7) solves the bilinear equa- 

tions with (ai = 0) 

A, = L - at, \i = L + en (i > 1), 
Ki = K — ai, Ki = K + ai (i > 1), 
a = K-L, (3 = K + L. 

Proof. First consider (J, i) with i > 2. Setting Zi = I2 = 1, Z = 0, x = 2;°, 63 = dj in 
(B.13), we have 

(6 2 -a i )F M ( a; o_ 1 + £ (6 2 - 1 ))F lil (a ; o+ £ (6r 1 )) 
- (61 -OiJFi,!^?.! +e(&r 1 ))*i,i (*° +*(&2 1 )) 
= (62 - &i)*i,i(*°-i)*i,i(*° + e(b^) + s{b^)), 

where the dependence on g e e D °° is suppressed. Due to (5.1) this yields (J,i) 
with 1 < J < 4 upon taking (h,b 2 ) = (L, K), (L, -K), (-L, -K), (-L, K). Next 
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we treat (J, 0) and (J, 1). Setting h = h = 1, 1 = 0, x = x^b^ 1 = in (B.12) =R 
(B.ll), we have 

b 2 F hl (x a 1 +e(b- 1 ))(F 0A T^F 0A . 1 )(x a 1 +e(b- 1 )) 
- hF^x* + e(b^))(F 0A T iFo,i;i)(*i + ^i" 1 )) 
= (& 2 - + + <b^))(F 0A t iFo^M). 

By setting (61,62) = (L,K), (L,—K), (—L,—K), (—L,K), and applying Lemma 
5.1, this becomes (J, l)r J and (J,0)t^ for 1 < J < 4. □ 

Step 2. We choose the elements g' £ e D '°° and g e e Doa related by (B.5) as 

N M 

(5.8) g' = cxp(j2 MfoMfc) + E ^(p'Ml'j)) , 

3=1 3=1 

N 



(5.9) g = exp(]T&,(<A (1) (f^ (1) *(-*) ~ V- (1) fe> (1) *(-ft)) 

M 

+ E^(^ (1) (^)^ 2) *(-^)-^ 2) (^)^ (1) *M))). 



On these elements we impose the reduction condition: 

(5.10) p 2 3 A(p,)A(-P3) = q*A( qj )A(- qj ), fiA^A^) = q') 

for all j, where the function A is defined by 



n-l 

a(p)=iiu-£) 

fc=2 fc 

in terms of 02, 03, . . . introduced in (5.4). 

Lemma 5.3. Under the condition (5.10), there exists h' G e D °° erne? /i S e- 000 £/iaf 
are related by h — t(h')Kt(h') as in (B.5) and satisfy 

Fi u i 2 -i{x + z n -i,y;w(g)) = F hM {x,y;h) 

for any x and odd y. 

Proof. Under the map u specified in (B.3), ip( a \p) and tp^*{q) are transformed 
into pip( a \p) and q^ 1 ip^*(q). Then under the time evolution Ade^"- 1 ' ), they 
are further changed into pA{p) Sal ip^ (p) and q~ 1 A(q)~ Sal ip^*(q), respectively. 
Therefore if g in (5.9) is denoted by e x , we get Ade^ 2 ™- 1 ,0 ) (cj(g)) = e x with 

(5.11) 

j J J J 
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Under the condition (5.10), we find 

X' = X\ b .^ v . tC ^ c '., 



p'M, 



In view of the definition (B.6) and H(x + z n -i,y) = H(x,y) + H(z n -i,0), the 
elements h! and h in question are obtained by replacing bj , Cj by bj , c'j in g' and g 
given in (5.8) and (5.9), respectively. □ 

Using bl and h in Lemma 5.3, we now modify (5.5)-(5.7) in Step 1 into the finite 
n version as 



(5.12) 

(5.13) r- 

(5.14) r$ 
Lemma 5.4. 



Xj = F hl (xf;g) \<j<n-2 (X = S,W,N, E), 



f 3 (x{;g') j — 0, 1, 

2<j<n-2, 
f n -j(x{;h') j = n-l,n, 

f (Fo,i + (-l)J*Fo,iii)(a:S;») 



(J= 1,2,3,4), 



^1,1(^-1; 5) 



(-1)" 



if L 



J =0,1, 
2 < j < n - 2, 
^^_ M;1 )(o:0_ 2 ; W (. 9 )) j = n-l,n. 



Fo,o(^_ i; — T~n—l^~n J =1,2, 3 ,4, 



T 1 T 2 
n n— 1 n'n— 



2 3 
T T 



I 3 2 
n 'n— 1 'n 'n— 1' 



r 4 r 3 
'n 'n — 1 



± T T 



3 T 4 
n ' n—1 5 



' n 'n—1 



n ' n— 1 ' 



2F ,o(a£_ 1 ;w(ff)) 
2iA'- 1 F_i,i i i(a:^_ 1 ;w( ff )) 

2^,0(^-1^(5)) 
2»L- 1 F_ lil;1 (xr_ 1 ;u;( 3 )) 

2F 0fi {x s n _ l -M9)) 
2iK- 1 F- 1 , 1;1 (xfl_ 1 ;cj{ S )) 

2F Qi0 (x%_ 1 Mg)) 
2iL- 1 F- lil;1 (x%_ 1 ;u(g)) 

Proof. The first relation follows from Lemma 5.3, Lemma B.l and (5.2). From 
Lemma 5.3, the left hand side of the second equation 2Fo t o(x^_ 1 ;u)(g)) is equal to 
2-Fb,o(zf; h). From (5.1), (5.3) and Lemma B.2, we get 

2F Q , (x?; h) = f (x{; h')h{x\; ti) + f {x\; ti)h(x\;ti), 
which is the right hand side. The other relations can be checked similarly. □ 

Theorem 5.5. The parameterization (5.12)-(5.14) solves the D„ bilinear equa- 
tions with (ai = 0) 

Aj = L - a,i, Xi = L + a,i (1 < i < n - 1), A„ = 1, 
Ki = K — a i} Ki = K + ai (1 < i < n — 1), n n = 1, 
a = K-L, [3 = K + L. 

Proof. Thanks to Proposition 5.2 we are left to show (J, n—1) and (J, n) only. Let 
us illustrate the proof for J = 1 since the equations for the other J follow from the 
same argument by taking the parameters b\ , 62 in the proof similarly to Proposition 
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5.2. In (B.13) set h = h = I = 0, {h,b 2 ) = (L, K), b 3 = a„_i, x = ar°_ x and 
g ^> u(g). The result reads 

K„-iF 0;0 (a;f_ 2 ; w(5))F 0)0 (a;^_ 1 ; - A n _iF , (x^_ 2 ; LJ (9)) F o,o{xf l _ 1 ; w(g)) 

Owing to (B.8), we know Fo i o(x^_ 2 ; u(g)) — Fi t \(x^_ 2 ; g) — and simi- 

larly F 0i0 (x%_ 2 ;uj(g)) = N n - 2 . Rewriting F 0j0 (a;i_i; w(g)), ^0,0(^-1; v(g)) and 
Fofl(x^_ 1 ;u>(g)) by Lemma 5.4, we find that the above equation coincides with the 
combination (1, n— l)r^_ 1 +(l, n)r^. The other combination (1, n— 1)t^_ 1 — (1, n)r r ^ 
can be shown by using (B.10) with Zi = Z 2 = / = similarly. □ 

5.2. Tropical vertex model. Consider the two dimensional square lattice £ 
equipped with the space-time coordinates (s,t) G Z 2 . The coordinate s (resp. t) 
increases rightward (resp. downward) and each vertical (resp. horizontal) line corre- 
sponds to s ^constant (resp. t =constant). On each edge of £ we assign an element 
in B so that those four surrounding a vertex obey the relation R(x,y) = (x',y'), 
where x, y, x', y' are the ones on the east, north, south and west edges. We call the 
resulting two dimensional system the tropical vertex model. Unlike the usual vertex 
models in statistical mechanics [B] , it is a deterministic system in the sense that all 
the edge variables are determined uniquely from their values on the northwest or 
southeast boundaries of £ . The tropical vertex model reduces to the D^-soliton 
cellular automaton [HKT1, HKT2, HKOTY] in the ultradiscrete limit, where the 
tropical R is replaced by the combinatorial R [HKOT] and B by the crystal 
[KKM]. From (4.4) it follows that all the elements in B on the same line possess 
the same level. We let the levels be l^ 1 (resp. k^ 1 ) on the ith horizontal (resp. sth 
vertical) line. 

The bilinearization of the tropical R attained in Section 4 leads to another for- 
mulation of the tropical vertex model in terms of tau functions. In view of Figure 
2 we duplicate each line in £ into a pair of parallel lines to form a new lattice £ 
Its unit structure looks as Figure 2 where we now replace {L,K) by (L(t),K(s)). 
To each face of C we assign a coordinate (s,t) which now takes values in (Z/2) 2 
so that r° : (s,t) G I?, N : (s,t- ±), W : {s - \, t), S : {s,t+ §), E : (s+\,t), 
t 1 : (s+ \,t- |), t 2 :{s- \,t- \), t 3 : (s - |,t+ \), r 4 : (s + \,t+ §). There 
are three types of faces depending on whether the coordinate (s, t) belongs to Z 2 , 
(Z + \) 2 or else. To the face of L at (s, t), we associate a tau function Ti(s, t) having 
the components 1 < i < n — 2ifs + teZ+i and < i < n otherwise. Then up to 

a boundary condition, the tropical vertex model is equivalent to imposing the 
bilinear equations on the tau functions around each face at (s, t) El?. 

Let us construct such system of tau functions {rj(s,i)} by making use of the 
result in Section 5.1, which may be regarded as a solution of our tropical vertex 
model. To each face of C at (s, t) e (Z/2) 2 we attach the time variables 

Xj(s, t) — x(s, t) + Zj, 

x( S , t) = v +j2 mt 1 )- 1 ) + E wr 1 ) + E £ wr 1 ) + E 

t'>t t'>t s'<s s'<s 
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where Zj is defined in (5.4) and r\ = i) is an arbitrary odd time. The sums extend 
over s', t' e Z. Let g,he e D °° , g' , ti G e D ~ be as in (5.12)-(5.14). Set 

T j (a,t)=F 1>1 (x j ;g), 1 < j < n - 2, s + teZ + 1/2, 
fj{xi\g') .7=0,1, 

Fi,i(^;<?) 2<i<n-2, («,t) G (Z+l/2) 2 , 
J n -j{xi;ti) j = n-l,n, 

'(Fo^ + i-iyiFo^ix^g) 3=0,1, 
Fi^(xj-i\g) 2<j<n-2, 
,( F o,o + ( ~ ) " 1 -F- 1 , l ; i ) (^-2 ; ^ (ff) ) 3 =n- l,n, 



= < 



= < 



(a,t) ez 2 , 



where 



irj(s,f), and the time variable y suppressed here is taken as explained 



in the paragraph preceding (5.5). Since x(s,t) obeys the recursion relation corre- 
sponding to (5.1), we conclude that {tj(s, t)} satisfies all the bilinear equations 
if the parameters X,K,a,/3 at (s,i) G Z 2 arc identified with those in Theorem 5.5 
with (L,K) replaced by (L(t),K(s)). In particular, the levels l^ 1 and kj 1 of the 
edge variables are determined by 



n-l 



n-l 



h = L{tf [] Wrf a?), k s = K{tf [] (K{sf af). 



Note that the formal sums over s', t' in x(s, t) make sense under a suitable choice 

Of 7]. 

Remark 5.6. Consider the composition of Figure 2 corresponding to RiR 2 Ri(x, y, z) ■■ 
{x',y',z') and R 2 R 1 R 2 {x,y, z) = (a;" , y B , «" ). 




According to the bilinearization, assign the common time variables and tau func- 
tions to the faces 1, 2, . . . , 7 of the two diagrams to represent the incoming state 
(x, y, z) € B x B x B. Then the time variables on the faces a, b, ■ ■ ■ ,e are also 
the same in the two diagrams. Therefore one has (x',y',z') = (x",y",z") for the 
associated outgoing state in B x B x B, which implies the Yang-Baxter equation 
R\R 2 R\ = R 2 R\R 2 . 
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6. Reduction to A < ^_ 1 and 

Here we introduce the tropical R for A^-i and Cn^ based on Remark 2.1 and 
explain how the results on the Dn^ case can be specialized to them. We note that 
these tropical R have more intrinsic characterization in terms of geometric crystals 
as argued in [KOTY] for D [ n ] . 

6.1. Tropical R for -A^n— l an< ^ ^n' '• ^ n tms sec tio n we denote the tropical i? 

for Z?! 1 - 1 by _R„ : B n x6„^ B„ x B n . Let B' = {a; = (xi, . . . , x n ,x n , . . . ,x"i)} be 
the set of variables and p : B' — > be the embedding defined by 

P((^l ) ■ • ■ 7 ^ni ^ni ■ ■ ■ i^l)) (Xl 5 • • • > X n , 1? X n j • ■ • j Xl) • 

The tropical R for A^-i is the birational map R' n : B' x B' — > B' x £?' defined by 
K(x,y) - (x',y') J R„ +1 (p(x),p(y)) = (p(x'), p(y'))- 

Let £?" = {a; = (xo,xi, . . . ,a; n ,5 n , . . . ,xi)} be the set of variables and p 1 : B" — > 
£>„+2 be the embedding defined by 

P ((^0; Xl , . . . , X n , X n , . . . , Xl)) = (xo, Xl , . . . , X n , 1, X n , . . . , Xl , Xo) • 

The tropical i? for c4 1} is the birational map R' r [ : B" x B" -» B" x B" defined by 

<(x,y) = (x',2/') ^i? n+2 (p'(x),p'(y)) = (p'(x'), pV))- 

The maps R' n and i?" are well defined due to Remark 2.1. They satisfy the inversion 
relation and the Yang-Baxter equation. We shall write R' n and R'^ simply as R' 
and R" when the rank needs not be specified. 

6.2. Bilinearization of R! and R" . First consider the tropical R' n _i for A^_ 3 . 
The bilinear equation for it is given by setting A n = K n = 1 and specializing the 
tau functions as 

(6.1) ri_ x =r J n < J < 4 

in the bilinear equations. The constraint (6.1) implies x n — y n — x' n = y' n = 1 
in (4.2)-(4.3), which is consistent with Remark 2.1 and Proposition 3.5. It makes the 
bilinear equations (J, n — 1) and (J, n) equivalent. We call the resulting subsystem 
{(J, i) I 1 < J < 4, < i < n — 1} of the case the A^l_ 3 bilinear equations. 

As for the tropical i?"_ 2 for C^]_ 2 i the bilinear equation is obtained by further 
setting Ai = Ai,Ki = ki and imposing the constraint 

(6.2) Tq == < J < 4 

on the A 2 ,i_ 3 bilinear equations. The constraint implies x\ = x\, yi = y 1 , x\ — 
x[, y[ — y[ in (4.2)-(4.3), which is again consistent with Remark 2.1 and Proposi- 
tion 3.5. It makes the equations ( J, 0) and (J, 1) equivalent. We call the resulting 
subsystem {(J,i) | 1 < J < 4, 1 < i < n - 1} of the A^_ 3 case the C^} 2 bilin- 
ear equations. Under these specializations Theorem 4.4 is still valid for R' n _ 1 and 

n n-2- 
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6.3. Solutions of bilinear equations for R' and R". Results in Section 5.1 
can be specialized to make (6.1) and (6.2) hold. First we deal with (6.1) for R' . 
Let us restrict the elements (5.8) and (5.9) to the case M = 1 and take the limit 

(6.3) q[ — > 0, p\ — > a m 

for some 2 < m < n — 1 keeping (l—Pi/a m )/q[ to be an arbitrary nonzero constant. 
It makes the latter of (5.10) void, leaving a\, . . . , a n -i as free parameters. All the 
tau functions in (5.12)-(5.14) are well defined in the limit as argued in the end of 
Appendix B.6. Moreover in (5.11), the coefficients p' 1 A(p' 1 ) / q[ and q' 1 /(p' 1 A(—p' 1 )) 
in the second sum tend to zero. Thus the element h! e e D ~ in Lemma 5.3 actually 
belongs to e B °°. Then the property (B.15) ensures that t^_ 1 = t;[ for 1 < J < 4 in 
(5.13) in agreement with (6.1). The remaining relation r°_ 1 = t" in (5.14) follows 
from (B.16). 

Let us proceed to R" , where the further condition (6.2) should be satisfied. The 
reduction from case to C^l 2 ^ s done by setting Vcj = in (5.8)-(5.9). Namely 
we restrict g' e e D °° and g e e D °° to g' € e B °° and g e e B °° , which reduces 
the relevant fermion theory essentially to the single component one. The former 
condition in (5.10) is still present although the latter becomes void. Then obviously 
h! e e B °° holds in Lemma 5.3, therefore (6.1) and (6.2) are valid due to (B.15), 
(B.9) and (B.8). 

Appendix A. Ui,Vi in terms of tau functions 

Let us prove Lemma 4.3. Our task is to substitute (4.2) into the functions Ui, Vi 
and simplify the result by means of the bilinear equations on the tau functions. 

Lemma A.l. 

11 /Jt? P±tlr ^ W-i „ 



Vn 1 P^n-l T n T , 



X n -1 + V n -l K n N n ^ 2 W n - 2 T^_ 1 T^' 

Proof. Apply (2,1),..., (2, rc-1). □ 

It is straightforward to check 
Lemma A. 2. For 1 < i < n — 1, t/ie ra£io (j/i • • • yi)/(xi ■ ■ ■ x^i) is equal to 

#(-D, Al ffff a (i = 2), Al -"^- 1 ^ (3<i<n-2), 
Ai---A„_ 2 W k - 2 (t 2 ) 2 t^_ 1 

132 2 ^ /' 

From Proposition 4.2, taking er* in Lemma A. 2 implies 
Lemma A. 3. For 1 < i < n — 2, t/ie ra£io (^i ■ • ■Xi)/{y 1 ■ ■ -I7j_i) is equal to 

m«=i), _ * iNi n 2 = 2), (3 <,< n - 2) . 

AiWi AiA 2 W2T 1 r 3 r 1 2 Ai • • • W^r 3 ^ " " 

In addition one has 



2 x 2 
n-1 T n-\ 



i-l'n 



Vl ■ Vn-lVn Al • • • A^-lTo^rX^T,^! 
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Proof of Lemma 4-3. First we consider Vo in (2.6). In view of (2.1) and (4.4), it is 
expressed as 

V = l~ 1 v+ + k^ 1 v^, 

n-2 

E » - , v 

xi---x n -i ( y n , 1 



2/1 ' 


■■Vi 






xi • • 


■Xi-l 


K.xi 


f v l 


• 


' ' X{ 






fi ' ' 


■Vi-i 







E yi ■ ■ ■ vi i j_ + j_ \ + yi • • - s/n-i / 

xi ■ • • T; i 1 X; 77 / m • • • nr~ o V ;r_ 1 

i=l 
n-2 

E 



2/1 ■ ■ ■ y n -iyn \x n -i y n -i 

Using Lemmas A.l and A. 2 we get 
To2o_ v = Ai AiTgr^ 1 ^ 2 y- 2 Ai ■ ■ ■ A^nfr/ Ai ■ ■ ■ A»_ 2 T 2 r°_ 1 r^ 

Rewrite t^t/ (1 < i < n — 2) and t^^t* appearing here by means of the bilinear 
equations (1, 1) — (l,n— 1). Then all but two terms cancel out, leading to 

fop ( 4 A i ■ ■ • Anr 2 r^ 

Similarly one can derive 

P T ( 4 , «1 ' ■•Kn-lToT^ 



W_ = 



Ai • • • A„_ir2 

from Lemmas A.l, A. 3 and the bilinear equations (3, 1) — (3, n — 1). Substituting 
these expressions of w± into (A.l), we obtain Vo = /3(fc — T o r o / '(^ aT o T o) as 
desired. Starting from this result on Vo, one can apply Lemma A.l and the bilinear 
equations (3, 1) — (3, n— 1) to the recursion (2.7) successively to derive the formulas 
for Vi, . . . , V„_i. Then the formulas for V^ 1 and V^ 7 * , . . . , V^ follow from (3.4)- 
(3.7) and Proposition 4.2. Finally substituting these results into (2.8) and using 
(4, 2) — (4, n — 2), one arrives at the desired expression for U\, . . . , U n -\- O 

Appendix B. Tau functions and free fermions 

Here we briefly recall the free fermion approach to the theory of tau functions 
based on [JM]. Lemma B.l and Lemma B.2 play an important role in the main 
text. 

B.l. Single component fermions. Single component charged fermions and the 
Fock states are given by 

bPj ,<Pk]+ = s ik, bPj , V'fe] + = bPj ,tPt]+ = ° U, k e z), 
Hp) = E ^ ' ^ ( p ) * = E 1 ' 

3 3 

(l\ = (vac|*?, 10 = *i|vac), 

' i>i-i ■ ■ ■ i>ii> , i>o (i>oi>* ■■■i>*-i, l >o 

= < i / = o , *r = { i i = o, 

k Vf •••V'-2^-i ^<o [v-iV>-2---^ ^<o 

ipj\vac) = 0, (vaclV'j = (j < 0), Vj|vac) = 0, (vac|V>j = (j > 0). 
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The neutral fermions are simple combinations thereof: 

+ i = T2 ' h = l 71 ' 

[<fij,(t)k]+ = m,<i>k]+ = (-!)%,-*;, [<f>j,4>k}+ = o, 

3 3 

(B.l) ^-10=^10=0 (Z = 0,1, j<0), (l\cf> j = (l\4> j =0 (Z = 0,1, j > 0). 
We prepare time variables and their functions: 

x = (xi,x 2 ,x 3 , ...), x = (xi,-x 2 ,x 3 , -x 4 , ■ ■ .), 

(B.2) e(a) = (a,££...) = -i(-a). 

We say x is odd if x = x. With the Hamiltonians 

m>l j£Z 

H'(X, y) = - Yl E {-l) m+1 {xi(t> m (t>-m-l + yi<f>m4>-m-l), 
1=1,3,5,... meZ 

the fermions exhibit the time evolutions 

Ade ff(a; V(p) = e« (x ' P V(p), kde H{x H*{p) = e -« (a: ^V*(p), 
Ade H '( x ' v U(p) = e^ x ' p U(p), Ade H ^ x '^4>{p) = e ^ v ^$(p) for (x,y) = {x,y), 
where Adg(-) = • 

B.2. 2 component fermions. 

= <W<^, [V^Uf ]+ = [^ a> ,4 0) *h = = 1,2, j, fc G Z), 

v (q) ( P ) = £ V, = 5>S Q) V J ', 

J 3 



" J V2 ' " J V2 

h(x, y) = + wvMSn, 

i>i j 

^ (a) |vac) = 0, (vacIV^* = (j < 0), ^ (a) *|vac) = 0, (vac|^ a) = (j > 0), 
(h,l 2 \ = <vac|v^**jf |/ 2)il> = *|>W|vac), 

where *| a) and * ; (q) * stand for and ^ with Vj.V'j replaced by ip ( "\ ip^* . 
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B.3. Algebras A^,, B' x , B x , D' x and D^. The algebras B'^ C D 1 ^ are de- 
fined within the single component theory as B'^ = {J2 a jk ■ 4>j4>k ■ +d \ B N 7 ajk = 
if \j + k\> N} and D'^ = {J2 a jk : <j>j4> k : + J2 b jk ■ 4>j<t>k ■ + E c jk ■ <t>j<t>k ■ +d \ 
B N, ajk = bjk = cjk = if \j + k\ > N}. The algebra Aoa in the 2 component theory 
is defined as = a ai , Pj : Vf^f h ■ +d | 3 iV, a aitfJj = if \i - j\ > N}. The 
automorphism a, k, uj of Aoo and the homomorphism l : D'^ — > Aoo are given by 

, (a) ■ i (a) i (a)* ■ , (a)* 

K : tpj I ^ llpj , Ipj I ^ —llpj , 

(B.3) u:^ a) ^^% ^-Vj-i*, 

The algebras B^ C Z?oo are defined as = {X e | <r(X) = X} and 

Boo = {-X' € Z?oo | 7r(X) = X}, where 7r is the projection to the first component, 
i.e., -K(ip^) — S a iipj a \ n(tpj a ^*) — S a iip^* . Thus B M actually stays within the 
single component theory. The map 

L + Kf.D^ —»£>«,, B^^Boo 

is an isomorphism and one has 

(B.4) L(H'(x,y)) + KL(H'(x,y)) = H(x,y) i£ (x,y) = (x,y). 

The set of well-defined group elements corresponding to these algebras will be 
denoted, by abuse of notation, by 

e D ~ : = { e Xl ■■■e Xh \ k > 0, X u . . . , X k e are locally nilpotcnt}, 

and similarly for A^, B^.B'^ and D'^. For g' e e D °° , we write 

(B.5) g = e^'H-PO = L ( g ') KL ( g >) g . 

B.4. Tau functions. We concern two kinds of tau functions: 

/l(x,»;fe) = <i|e H '( x,, '^|0 for(x,y) = (i,y), h e e D ™ , 
F hM {x,y;g) = {hM\e H ^g\h-l,h + l) g e e A °° . 

Fi lt i^-Q[x,y; g) will simply be denoted by Fi l j 2 (x,y; g). They enjoy the symmetry 

(B.7) F hM (x,y-g) = (-l^+^F^^.^y^ig)), 

(B-8) F hM {x,y-g) = F ll - 1M - 1 . t i{x,y;u{g)). 

Clearly we have 

(B.9) F hM (x,y;g) = if Z ^ and g E B M . 

Fi uh .i(x,y.g) will be denoted by F hM {x,y) or F hM {x) for simplicity. 
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B.5. Bilinear identities. For any g £ e Aa ° one has the hierarchy of equations 

/ ^(-y^k 11 ' 1 ' 1 ' 1 ^^'^^-^!^ - e(k- 1 ),y)F l , +1 ,,,,_ 1 ( X ' + e(k- 1 ),y') 

+ j ^ kh ~ l ' 2 ~ leiiy ~ y '' k)Fh ' l i- 1 * ix > V ~ ^-^F^+^ix', y' + eik- 1 )) = 0, 

where — > V — I > 1' 2 — h + 2 and the integration is taken along a small contour 
around k = oo. This is eq.(4.4) in [JM]. Let 61,62,63 be nonzero constants. We 
take y — y' and suppress the dependence on it. The above equations contain the 
following: 

(B.10) (62 1 - b^F^.^x + e(b^) + e(6i- 1 ))F il _ 1 , i2+1;(+1 ( a; + e(b^)) + eye = 0, 

(B.ll) (&2 1 - bs^F^ix + e(b^) + e{b^))F h _ 1M ,i +1 {x + e(b^)) + eye = 0, 

(B.12) (62 1 - b^F^aix + e(b^) + e (6J 1 ))F, 1 _ 1 , !2;; (^ + e^ 1 )) + eye = 0, 

(B.13) (62 - M^i.wOe + e^ 1 ) + e(&3 1 ))^i,! 2 ;i(^ + e^ 1 )) + eye = 0, 

where +cyc means the additional two terms obtained by the cyclic permutations of 
61,62,63. They also contain 

F .i(x + e(6 2 - 1 ))F .o(x + e(6r 1 ))- F ,i(a; + e(b^))F a . (x + e^ 1 )) 

( B - 14 ) 11 11 

= (6r 1 -6 2 - 1 )Fi A -i(x + e(6r 1 )+£(6 2 - 1 ))F_i 4; i(x). 

When both a; and y are odd, F; lj ; 2; ;(a;, y; 3) with 17 S e D °° can be expressed in terms 
of fi(x,y;g') as noted in eq.(7.6) in [JM]. We recall this fact in 

Lemma B.l. Let g £ e D °° and g' £ e D =° be related as in (B.5). Suppose that x 
andy are both odd. Denoting F hM (x,y; g) and fi(x,y;g') by F hM {x) and fi{x), 
one has 

F h i{x) =F 0fi {x) = f (x)f 1 (x), 
F , 1 (x)=F lfi (x) = ±(f (xf + f 1 (x) 2 ), 

Fo.iAx) = -Fi,Oi-i(aO = \{k{xf - fi(xf ). 

Proof. The first equality in each line is due to (B.7). Since <j>Q = = 1/2, we 
may write Me H ' {x ^g' = X + X X (j) Q + X 2 <t> +X 3 (f> a <j>o, where Xj = Xj(<j>,4>) 
contains <f>k,<f>k only for k ^ and its order is even (resp. odd) for j = 0, 3 (resp. 
j = 1,2). Due to (B.l) we find fi(x;g') = (X Q ) + {-l) l i{X 3 )/2 for I = 0, 1, 
where (Xj) := (0\Xj\0) = (1|X,-|1). On the other hand from (B.4) and (B.5) it 

follows that Ade H l x >rig = {X' Q + X^ + X 2 (f> (2) + X^Vo^X^o + X^ + 
0^0^ )> where X'- and X'- are obtained from Xj by substituting 
the 2 component fermions into Xj as X' 3 = Xj((f>^\ <^ 2 )), X] = X 3 {4> {1 \ 4> {2) ). 
Due to the property (B.l) for each component, we get F fi{x;g) — (0, 0|(XqXq + 

X^X^U^^^MO) = (X Q ) 2 + (X 3 ) 2 /A=((X )+i(X 3 )/2)((X )-i(X 3 )/2), 
proving the first relation. The other relations can be shown similarly. □ 

For odd x and y, it is clear from the above proof that 

(B.15) fo(x,y,g') = fi(x,y;g') if 5 'ee B ~. 

Lemma B.l admits a generalization to 
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Lemma B.2. Let g e e D °° and g' e e D ~ be as in (B.5). Suppose y is odd and 
x = x + e(c _1 ) - e(c _1 ). Denoting F h j 2[ i(x,y; g) and fi(x,y;g') by F h .i 2 ,i(x) and 
fi(x), one has 

-2iFo°viU } = MX + £{c ~ 1))h{x i{c ~ 1]) ±h{x + < c ~^ x ^ c ~^ 

Proof. In (B.10)-(B.12), set h = -b 2 = c and b^ 1 = 0. Writing x + = x + 

e(c _1 ), a;_ = x — efc" 1 ), one has 

F 1A {x-)F 0A (x+) + F 1A (x+)F 0A {x-) = 2F 1 . 1 (x)F . 1 (x), 
F 1A (x_)F 0A . 1 (x + ) + F ltl (x+)F 0>lil (x-) = 2F 1A (x)F 0A . 1 (x), 
F lfii _ 1 (x-)F , 1 (x + ) + Fi A -i(a;+)Fo,i( a; -) = 2F lfi ._ 1 (x)F 0tl (x). 

By applying (B.7) together with a(g) — g, the right hand sides here are equal to 
2^o, o(x)F 0t i(x), 2_F , o( x )Fo, i;i(%) an d —2F ,ia(x)F .i(x), respectively. Since both 
x + and are odd, the left hand sides are expressed in terms of fo(x±), fi(x±) by 
Lemma B.l. Solving the resulting three equations with respect to Fo t o(x), Fo,i(x) 
and Fo,v,i(x), we obtain the formulas in question up to an overall sign, which can 
be fixed by considering the smooth limit g — > 1. The formula for F-i.i-i(x) can be 
derived similarly by using (B.14). □ 

Note that in the limit c — > oo, the first three relations in Lemma B.2 reduce to 
Lemma B.l. 

B.6. Explicit formula. Let us write down the explicit formula for the vacuum 
expectation values (B.6). To simplify the result we introduce the parameters 
Pi, ■ ■ ■ ,P2N+m and qi, . . . ,qM to modify the parameterization of g' e e D ™ and 
g e e D °° in (5.8) as 

N M 



g' = cx p(X1 b i<i>(pi)<i>(pi) + cj<j>(pj)4>{qj)} , 
»=i j=i 

JV 

g = exp (J^ bi (Pi )4> W * (-Pi) - (~Pi )) 

i=l 
M 

+ ^ Cj (^( 1 )(p j )^ 2 )*(^)-^ 2) fe-)^ (1) *(-^))), 



where i = 2N + 1 - i and j = 2N + M + 1 - j. Given J C {1, . . . , N} and 
J C {1, . . . , M }, we shall write I = {i \ i £ 1} and J = {j | j e J}. The cardinality 
of a set 7 is denoted by |7|. For any sets K and K' we set 

Aiw= n (pm±p-). 

LL,V£K, (A<V 

**Up) - u Ak{p)Ak (p { ~? p ) n ?r n 

23 



and similarly for q. Let x and y be odd time variables. Then we have 

fi(x,y,g') = X) ei , J 2-l^-iJ| t/ ( g ) Cj(a;ii>) ^jP)Mg) , = 0,1, 



£(p)Aj(g) 

fl IJI PVI 

A _ = /U/UJ, ejj 



1 |J| even 

|J| odd, 

iei 

cj{x,y) = Y[cj exp(£(x,p~j) + £(y,qj)), 

where the sum ^ 7 j extends over all the subsets I C {1, . . . , A/"} and J C {1, . . . , M}. 
For any integers l\,h and Z, Fi lt i 2 - t i(x,y; g) with odd y (x is not restricted to be 
odd) reads 

F h,l2;l( x ,y,9) 

= E (-l) |7 ' l+(|7|+|/ ' l)|J| &/W^Wc7(^,y)c^( a; ,y)A^^(p)A^- i ( g ), 
i,i',j,j' 

K = IUl'Uj, K' = I'UIU J', 

iei' 

c '.i'( x ,y) = n c J ex p(-^( a; '-^) + 

where the sum ^ 7 7 , j J; extends over all the subsets 1,1' C {1, . . . , AT} and J, J' C 
{1,...,M} such that I J' I - |J| = I. 

Consider the dependence of fi{x,y;g') (I — 0,1) and Fi lt i 2; i(x,y; g) on q\ when 
M = 1. Apart from the factors c.j(x,y) and c'j,(x,y), ft is independent of gi and 
Fh,l 2 ;l is proportional to the single power (fr 2- . Therefore all the functions ft and 
Fli,h;l appearing in (5.12)-(5.14) are well defined in the limit q\ — > with M = 1. 
In particular from (B.8) one has 

(B.16) limF_ 11 .- L (x,y;ui(g))= lim F 0>2 ;i(x,y;g) = 0. 

gi ->0 gi~»0 

These facts are used in Section 6.3. 
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